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Instanton sector of correlated electron systems as the origin of populated pseudo-gap 

and flat "band" behavior: analytic solution 

S.I. Mukhiij:] 

Theoretical Physics Department, Moscow Institute for Steel & Alloys, Moscow, Russia 

(Dated: July 26, 2008) 

Finite temperature instantons between meta-stable vacua of correlated electronic system are 
solved analytically for quasi one-dimensional Hubbard model. The instantons produce dynamic 
symmetry breaking and connect metallic state with the dual vacua: superconducting (SC) and 
spin-density wave (SDW) states. The instantons spread along the Matsubara's imaginary time 
and possess the structure similar to the coordinate-space solitonic lattices previously discovered in 
quasi one-dimensional Peierls model. On the microscopic level the inter-vacua excursion is described 
by mutual transformations between the "resonating quartets" of the couples of electron-hole and 
Cooper pairs. Spectral properties of the electrons in the "instantonic crystal" reveal pseudo-gap 
(PG) behavior, with finite fermionic density of states in the center of the PG and "flat-band" outside 
of it. Analytically derived inverse temperature scaling of the pseudo-gap and the densities of the SC 
and SDW condensates is discussed in the context of ARPES and STM data in high-Tc cuprates. 



It is for long time now that competing orders are 
considered as one of the core properties of underdoped 
high-Tc cuprates [l|,[2l- Nevertheless, the nature of the 
non-superconducting order remains a mystery. Espe- 
cially challenging is the pseudo-gap (PG) phenomenon 
observed in high-Tc cuprates [3[. Here we study a novel 
possibility, which exploits a well known concept of high- 
energy physics, an instanton, that connects vacua of the 
many-body system possessing different symmetries |j]. 
The instanton sector of the Hilbert space of the many- 
body system is known to possess dynamically generated 
masses of the excitations in the system [J| . We find that 
this mechanism in the correlated electron system has in- 
teresting peculiarity: the "mass" is acquired by fermionic 
Landau quasi-particles inside the pseudo-gap, but neither 
right in its center, nor far away from the PG region of 
the spectrum. In the center of the pseudo-gap the quasi- 
particle dispersion remains linear in momentum ("rela- 
tivistic"), but acquires just strong enhancement of the 
Fermi velocity. Outside the PG region the spectrum re- 
mains not renormalized. As a result, in the cross-over 
region a "flat band" forms in the fermionic excitation 
spectrum, see Figs. 1,2, with a weaker dependence on 
momentum. We demonstrate these properties using a 
(quasi) one-dimensional model of electron system, which 
is especially convenient tool for understanding the major 
effects that occur, since it proves to be analytically solv- 
able and, hence, provides in-depths information. We also 
discuss predicted measurable characteristics of the elec- 
tronic system that undergoes symmetry excursions into 
the (dual) meta-stable vacua described respectively with 
superconducting (SC) and spin-density wave (SDW) or- 
ders. These vacua are chosen as the most relevant ones 
since ARPES, STM and INS experiments [1, H, @, 0, 3 
reveal PG behavior in the anti-nodal region of the Fermi 
surface of doped high-Tc cuprates with coexisting mag- 
netic and superconducting fluctuations well above the 



superconducting transition temperature even in the opti- 
mally doped compounds. From our present point of view, 
supporting evidence for the discussed instanton scenario 
also comes from the spectral structure of the magnetic 
fluctuations measured in doped high-Tc cuprates [9| near 
the anti-ferromagnetic points of the Brillouin zone. The 
structure reveals "rod-like" behavior of the intensity as 
function of frequency. This, when Fourier-transformed 
back to time, points to strong unharmonicity of the mag- 
netic fluctuations, thus reminding an "instanton hair- 
comb" rather than commonly considered RPA-like sinu- 
soidal fluctuations around an equilibrium. A complete 
theory should describe transitions between different sym- 
metries representing the 5*0(5) group: charge- spin- and 
superconducting orders, relevant to correlated electronic 
system. Nevertheless, we believe that SC and SDW or- 
ders could be less coupled to the lattice degrees of free- 
dom (due to charge transfer in cuprates) than CDW or- 
der, thus forming most probable end-points for the "fast" 
instanton routes, of relevance for the 10 -^ 100 meV scale, 
where the "rod-like" spectral intensity of the magnetic 
fluctuations is observed [9[. As we address behavior at 
flnite temperatures, an instanton propagates in the Mat- 
subara's imaginary time [lO| and is contained in the parti- 
tion function Z of the quasi one-dimensional many-body 
system [i,[ll|: 

Z ^ I V^V^VMVAexp{~{SF + Sep)} (1) 

f/3=l/T 
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Here we put the Boltzman's and Plank's constants to 
be unity (fc^ = h ~ 1). Hamiltonian of the interacting 
electrons in the decoupling Hubbard-Stratanovich (HS) 



fields M{t,x) and A{t,x) has the form (compare [2]): 

HeF = ^a i -i-^ ] <7z'^a + M{t, x)^aa+'^a + 

M{t, x)^a^-'i'a + A(T,x)(-\i' + J§_,x + ^-J^ + ^x) + 

A(r,a:;)(-1'_.x*+,t + *+,X*_,l). (4) 

here spin indices are either indicated by up and down ar- 
rows or in the summs by index a. The PauU's matrices a 
are defined on the pseudo-spin 1/2 space with -1-1/2 and 
— 1/2 projections distinguishing between right and left 
moving electrons in the one-dimensional system with co- 
ordinate X. The operators ^ and VP create and annihilate 
fermions and are taken in the Heisenberg's representa- 
tion with the Matsubara's imaginary "time" t. Partition 
function in Eq. ([1]) contains only (fermi)bosonic fields 
that obey (anti)periodic conditions, with the tempera- 
ture T defining the period /3 = T-'^ [uj: 

^(r-f-/3),*(r + /?) = -*(r),-*(r); (5) 

M(T + /3),A(T + /3) = A/(r),A(T). (6) 

Integration over fermionic (Grassmannian) fields ^ and 
^ in 111) could be (formally) done exactly, provided that 
one knows the full spectrum {«„} of the quasi-energies 
of the related Floquet equation: 

dr^n + H,f{M{t), A(r))$„ = (a„ + iTip-^)<^n (7) 

that follows from the Floquet 's theorem for the fermionic 
wave-vunctions in the "time"-periodoc bosonic fields 
M(r),A(r) (analo g of the Bloch's theorem for space- 
periodic potentials) [ll|: 



^/^*„(t) 



^$„(r); $„(r + /?)-$„ (r) (8) 



Since the indices a„ are known, partition function in Eq. 
(dl) can be expressed as: 



Z = I VMVAexp{-SF} Y[ c/^ 



(9) 



In the approximation studied below, the path-integration 
in Eq. ([9|) is substituted by the saddle-point value in the 
Hilbert space of the HS fields M, A, that is found from 
zero variational derivatives equations: 



'JM{T,x),A{r.x)- 



Sf — y^ In 



H^)}\-om 



While obtaining the self-consistency equations in the ex- 
plicit form below we were guided by the following math- 
ematical and physical facts and considerations. Firstly, 
the functional derivatives of the eigen- values q;„ with re- 
spect to the potentials M{t,x) and A{t,x) in pop are 
readily calculated using first-order perturbation theory 
applied to Eq. dH): 



Secondly,the Hamiltonian matrix Jfei?(M, A), as it fol- 
lows from Eq. ([4]), is a 4 x 4 matrix defined on the 
quartets of functions u^,u^,v^,v_(t,x). The latter 
are considered as the amplitudes of the right- and left- 
moving electrons and holes coupled by SDW and SC 
fields, which are searched for in the form: A/(r, x) = 
M{T)exp{ikFx} + M*(T)exp{—ikFx} and A(t), kp be- 
ing the one-dimensional Fermi-momentum. It proves to 
be, that in the static (i.e. r-independent) case, for any ar- 
bitrary modufi ratio of the SDW and SC orders: |M|/|A|, 
and for arbitrary phases 4>w and 0s of these order pa- 
rameters, the matrix HeF{M,A) can be diagonolized, 
provided that the u, v amplitudes obey the following con- 
ditions [21]: 



v± = e*(=^'^»--*^)Mzp; v± 



J(±4>\v-<Ps)^ 



(12) 



Then, it is logical to assume that the instanton passes 
between SC (|M| = 0) and SDW (|A| = 0) vacua 
along the trajectory in the Hilbert space curbed by the 
{u-i_, u_, w-i_, w_}-quartets, that obey conditions ()12p at 
any moment inside the imaginary time period < r < /?. 
Under these conditions the Hamiltonian 4x4 matrix is 
2 X 2-block diagonal, and equations defining the eigen- 
functions ^n{T) in Eq. ^ for each 2x2 block read: 

{a^dr - A± {t)cf+ + A± {t)cf^ + vkl) Uk^O (13) 

where {uk] ~ {u+,u^) are Fourier-components of 
the slowly varying in space (x-coordinate) amplitudes 
u± (r, x) of the right/left-moving fermions with the wave- 
functions: u = u^{T,x)exp{ikFx) + u^{T,x)exp{—ikFx); 
V- is Fermi velocity. The combined HS fields A± are 
defined as follows: 



A±(r) = (|M(r)|±|A(T)|)e» 



(14) 



The choice of the sign + or — inside A± distinguishes 
the two different 2 x 2-block parts of the subdiagonalized 
complete 4x4 Hamiltonian matrix. Below we assume 
A± to be real. Then, a unitary transformation: f± = 
(2)"i/2(u+ ± u_) transforms Eq. ^ into: 



{dr + vk&a; - A±(r)(7^) fk = 0; 



(15) 



where {fk}'^ ~ (/+, ./-). In the case of the self-consistent 
instanton solutions with real fields A± = A± found be- 
low, the (Matsubara's)time-conjugated functions obey 
the following symmetry relations: 



f±{k)^hi-k), 



(16) 



d 



M,A{T,x}<^n 



a„ = $„(r,x){aA/,Ai^eF(Af, A)}$„(t,2;) (11) 



that can be checked directly. Functions f± satisfy equa- 
tions following from Eq. ^T5\\ : 



(17) 
(18) 



id\-Q±{T)~{vkf)f±=0; 



and equations for the conjugated functions f± can be 
obtained from Eq. (J17p by the siniuhaneous interchange 



/(fc)d 



,f{—k) and Q±{t) -^ Q^{t). It is obvious 



from Eqs. pT]) . (fl5|) that the y d iffer from the correspond- 
ing equations for the Peierls [1^, ll3| and one-dimensional 
Hubbard model [lj| only by the opposite sign in front of 
the square of the bare dispersion vk and by the inter- 
change of the space variable with the imaginary Matsub- 
ara's time r. This latter interchange defines the period 
of the A±-functions to be equal to the inverse temper- 
ature (5 — l/T, while in the Peierls model a real-space 
period is proportional to 27r/e, i.e. to the doping concen- 
tration e away from the half-filling of the bare electronic 
band. Hence, a decrease of the temperature corresponds 
to "underdoping" . 

An explicit expression for the functional derivative of the 
Floquet indices a„ with respect to combined HS fields 
A± is readily derived from Eq. (fTS]) by means of the 
first-order perturbation theory: 



5A±(r)afe 



{f+{k)f^{-k)- 



f-{k)U{-k)}t 



(19) 



where the last equality is obtained using the symmetry 
relations in Eq. p6)) . and superscript ^ distinguishes 
solutions of Eq. p?)) with the two different HS fields 
A± defined in (fT4|) . Equation (fT9|) can be now substi- 
tuted into equations (fTO|) . resulting in the following sys- 
tem of self-consistency equations for the two instanton 
fields A± (linearly related with the SC and SDW fields: 
A = 1/2(A+ - A„) and M = 1/2(A+ + A_)): 



^th 



ak~ 



{/+(fc)/_(-fc)-/_(fc)/+(-fc)}± = 



A+(r) f[/-i + W) ± A_(r) (u'^ W^m 



We found that Eqs. HZ]), HI]) and ^ possess self- 
consistent instanton solution that can be obtained di- 
rectly from the self-consistent soliton solution found pre- 
viously [1^, ll3[. Thus, the former solution can be ex- 
pressed via the Weierstrass's V and Jacobi's snoidal 
A2sn(Ai;fci) elliptic functions, that also described soli- 
tons in the one-dimensional chains. Only now the 
Matsubara's imaginary time plays the role of space- 
coordinate, hence, leading to the following self-consistent 
instanton solution of the equations (fT7|) . (|18p and ([20|) 
(compare 12|, [iSj): 



/±(T,fc) = ^exp|TC(fc)^'0±(t±)dr'|(21) 



w±{t) = v/e-t±; (j){t) = 



1 



(22) 
(23) 



e — t± 

t±^'P{T + T±); T±^UJ2.3+To; 
3 

C{k) = l[^{e-e,); e ^ ei + {vkf (24) 



Here TV = e — "P is a normalization constant (constant V 
is explicitly expressed below), ei > 62 > 63 are the three 
values of the Weierstrass function V at its three half- 
periods wi, LO2 and ^3 ~ — (wi + LJ2), that are related 
with the parameters Ai^2 and ki of the corresponding 
Jacobi function A2sn(AiT; fci). The latter describes the 
instantons A±(r): 



A± = A2 sn(Afr;fcj^ 



{A±}2 + {A±} 



±12 



'2,3 



2 2 ' 1 



A| 
A? 



(25) 



It is well known that both elliptic functions V and A± = 
A2sn(AiT;fci) share the same period 2uJi as functions 
of the real variable r, provided they are related in the 
following way [15[: 



2V{t + T±) + ei = A^±) ± a,A(±) 



(26) 



Hence, the Matsubara's time periodicity conditions ^ 
lead to the following self-consistency equation (two equa- 
tions appear if we retain the indices distinguishing be- 
tween the two instanton fields Aj-): 



1 „ _ 2X(fc) 

^ = ^ = ^^1 = 7 u72 

T (61-63)^2 



;fc^ = 



62 -63 
61-63' 



(27) 



where T is the temperature, and K{k) is the complete 
elliptic integral of the first kind. The remaining two self- 
consistency equations follow then from Eq. (jSO]) . that 
due to Eq. ([21]) become remarkably simple: 



2A±(T)^t/i 



f Oik'- 



\ 2 



vk 



V^ 



^*'^'[v^¥l}^^-'^'ih-Wl 



(28) 



where V is half-period averaged Weierstrass function: 



V = ax = 61 — (ei 

Jo ^1 



^'>f|'-' 



where i*^(fc) is the complete elliptic integral of the second 
kind. It is possible now, using definition of the Floquet's 
indices from Eq. ([5]) and solutions for the wave-functions 
from Eqs. (PT|) - (|M|) . to find an explicit analytical ex- 
pression for the spectrum of fermion quasi-energies in 
the presence of "traversing" instantons: 



at = Cik) 



P=2uJi 



2V{i 



-dT 



(30) 



i=l 



The well known properties of the Weierstrass function 
[15j and of the complete elliptic integral of the third kind 



n(n, k) enable us to find the analytic answer in a compact 
form (omitting the superscript for the sake of clarity) : 



ak = 2e/t 



1-e 



1/2 



n 



:,k\ (31) 



vk 



Cfe = 



2TK{k) 



(32) 



This expression is remarkable, since it reveals the pseudo- 
gap with "localized energy" states inside of it, and a 
"flat-band" behavior between the gap edge and the "free- 
fermion" spectrum, see Figs. 1,2. These features re- 
mind well known data measured e.g. by ARPES in high- 
Tf, cuprates. The corresponding asymptotic expressions 
read: 



Oik 



vk 0(1) 



au ~ 2ln 



T k'K{k) 
( ATK{k) 



'—] vk< k'TK{k) 
vk 



\ vkk' 



k' < 



TK{k) 



-^ < 1 



Qfc w — ; TK{k) <vk< ep- 



(33) 
(34) 
(35) 



where ep is of the order of the bare electron energy cut- 
off, and k'^ = 1 — k^, to be found self-consistently below. 
Using definition of e from ([M)) and equation (|29p we 
rewrite self-consistency equations (P5| in the form: 



(36) 



A±(T)7±=A+(T)(5+±A_(r)<5_ 
where we introduced notation: 

S+^iU-'±\g\-'); 



and the essential coefficients "f± are: 
' 12/ r.,1.^2 ^/i/T 

k 



7± 



2j2th 



2 J {vk)^+4{T)^K{k±)E{k±) 



(37) 



(38) 



We still have choice with respect to the "initial times" 
Tq in Eq. §^ since A±(t) = Afs7i{Af{T + T^);kf). 
Since the self-consistency equations Eq. ((36|) fulfill in 
each moment of "time" r there are only two irreducible 
possibilities: Tq = and Tq — Tq ~ wi/Ai. In the 
first case A±(t) cx sn{Af{T);kf), and in the second 
case A±(t) cx ±sn(Af (t); fcj^). The second case proved 
to have no solution with the pseudo-gap property of the 
fermionic spectrum. Hence, we consider below the first 
case, i.e. Tq = 0. It follows directly from Eqs. (|33p - ((35|) 
that essential (logarithmic) contribution to the sum in 
Eq. (1551) comes from the momenta interval: k'TK{k) < 
vk < ep- This leads to the following result: 



^^ ~ —In ( ep/2T\lK{k±)E{k^) 



(39) 




20 30 
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FIG. 1: Pseudo-energies spectrum for different values of the 
coupling strength parameter k' as given in Eq. (|44p . see text. 



Equations (P5)) - (P5)) indicate that the pseudo-gap behav- 
ior is present provided that k' « 1. The self-consistent 
solution possessing this property reads: 



A(^) = ^^Mir) 



(40) 



M{t) = Am ( 1 - ^^^77^ ) sn{AK{k)TT- k); (41) 



2U 



^-^i^-{-¥f^'^<^* 



(42) 



k' = 4exp{-Qy, e = (g) exp |-^ \ (43) 



Am = 0, T > T*; T* cx epexp { -^ 

m 



(44) 



Substituting k' from Eq. (gH) into Eqs. (I231)-(13S]) we 
evaluate the pseudo-gap ApQ-. ApQ cx Am ^ 1/T. 
In summary, we see from Eqs. (pn i) - (|44p that when elec- 
tron coupling constant \g\ is close to on-site repulsion 
U the co-existing superconducting and spin-density wave 
fluctuating condensates (instantons) form below the tem- 
perature T* that depends on superconducting coupling 
l^l, with both condensates densities being inverse pro- 
portional to temperature. The density of the supercon- 
ducting condensate A vanishes at the point \g\ = U. 
Also, a pseudo-gap in the fermionic energy spectrum ex- 
ists with the magnitude Apc proportional to the den- 
sity of the magnetic condensate Am- In the center of 
the pseudo-gap there are fermionic states with the den- 
sity of states that decreases exponentially with temper- 
ature according to "quadratic Arrenius" law ^ k' ~ 
exp{~'Const ■ {£p/2T)'^}, where const = exp{—2£p/\g\. 
These in-gap states are also observed experimentally 3| 
in high-Tc cuprates. 

Also in qualitative accord with the ARPES experiments 
in these compounds [3|, |5|, |6[ we find at the edge of the 
pseudo-gap a "flat band" in the fermionic dispersion, see 
Fig. 2. 

As it follows from Eq. (|33p . when substituting into it Eq. 
(im, the density of states in the center of the PG grows 
with decreasing coupling strength Ig] ~ L'^ (and thus in- 
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FIG. 2: Density of states (DOS) for the same values of k' as 
in Fig. 1. 



creasing parameter k' in Eq. (011)), while the PG am- 
plitude decreases simultaneously. When combined with 
Hall effect data in cuprates discussed in [16|, the in-gap 
states and PG behavior found in the present work may 
be tolerated qualitatively with the PG behavior inferred 
from Hall experiments if we assume that couplings \g\,U 
decrease with the hole doping of the high-Tc cuprates. 
Finally, we speculate, that if the instanton contribution 
to the partition path-integral Eq. ([T]) becomes negligable 
with respect to the contribution of the static SC or SDW 
order when the temperature T lowers, the ground state 
is either superconducting or antiferromagnetic. On the 
other hand, another scenario to be considered is when the 
system freezes while remaining in the instanton sector of 
its Hilbcrt space. The set of parameters that enables such 
behavior will be investigated in the forthcoming work. 
The author appreciates stimulating discussions with 
the colleagues at the Lorentz Institute for Theoretical 
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